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A model independent determination of \V u b\ using the global q 2 dependence of the 
dispersive bounds on the B — > irlv form factors 
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We propose a method to determine the CKM matrix element \V u b\ using the global q 2 dependence 
of the dispersive bound on the form factors for B — > nlv decay. Since the lattice calculation of the 
B — > -kIv form factor is limited to the large q 2 regime, only the experimental data in a limited 
kinematic range can be used in a conventional method. In our new method which exploits the 
statistical distributions of the dispersive bound proposed by Lellouch, we can utilize the information 
of the global q 2 dependence for all kinematic range. As a feasibility study we determine \V u b\ by 
combining the form factors from quenched lattice QCD, the dispersive bounds, and the experimental 
data by CLEO. We show that the accuracy of \V u b\ can be improved by our method. 

PACS numbers: 12.38.Gc (temporary) 



I. INTRODUCTION 

Precise determination of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements from the B, D, and 
K decays is one of the major goals in flavor physics. By 
measuring both the sides and the angles of the unitar- 
ity triangle, one can test the consistency of the standard 
model, which can either verify the standard model or 
probe a signal of new physics. Although useful for the 
consistency check against sin(20i), the CKM matrix ele- 
ment \V u b\ is one of the most poorly known quantities at 
present. It would be important to reduce the uncertain- 
ties further. 

There are two ways to determine | V u b , i.e. the determi- 
nation from the inclusive semileptonic decay B — > X u lv 
, and the determination from the exclusive semileptonic 
decay B — * -kIv or B — > plv. Since both methods suffers 
from different systematic errors from the experiment and 
the theory, having independent results from the inclusive 
and the exclusive processes are necessary for the reliable 
determination of \V u b\- 

In the former method, the hadronic matrix element in 
the inclusive process is the operator product expansion 
(OPE) in l/m&. It is suggested by Bauer et al. [l| that 
with an appropriate choice of the combined kinematical 
cut (m 2 x and q 2 ) to remove the charm background, the- 
oretical errors from various higher order corrections in 
OPE are controlled so that \V u b\ can be determined at 
the level of around 10% or below. 

The latter method requires the computations of the 
form factors which describe the hadronic weak matrix 
elements for the exclusive processes. Lattice QCD is a 
promising tool for this purpose. However at large recoil, 
the cutoff effects of order O(aE) becomes non negligible 
where a is the lattice spacing and E is the recoil energy of 
the daughter mesons such as 7r or p. Therefore the precise 
computation of the form factor is limited only for large q 2 
region in the lattice QCD method. Although B factory 
experiments can measure the differential decay rate at 
high q 2 regime, they are statistically limited while the 



rich experimental data for the low q 2 region will remain 
untouched. It would be ideal if one could make a precise 
prediction of the semileptonic decay form factors for the 
whole q 2 range. 



Lellouch proposed a statistical method in which one 
uses the form factor values with theoretical errors for 
B — > irlv at high q 2 from lattice QCD to give a distribu- 
tion of the dispersive bound for the form factors for the 
small q 2 region 2] . Although the idea is attractive, with 
his original idea one could only obtain a loose bound for 
the form factor. Later several authors pointed out that 
this bound can be improved with additional inputs near 
Qmax M ■ Recently CLEO gave the q 2 dependence of 
the B — > ttIv decay Although the overall normal- 
ization for the form factor is unknown up to \V u b\, the 
CLEO results give a strong constraint on the q 2 depen- 
dence of the B — * nlv form factor. In this paper, we 
propose to use the information on the q 2 dependence of 
the B — > irlv from the experiment to restrict the statis- 
tical distribution of the dispersive bound. We show that 
this additional input leads to a better determination of 
\V ub \. 



This paper is organized as follows. In [n] we explain 
our basic idea for \V u b\ for determination using lattice 
data, dispersive bounds, and the experimental data. The 
explicit method how to combine those data explained in 
detail in section IIIII In section IIVI our results on the 
improved bound for the form factor as well as \V u b\ are 
presented. In section^ we discuss the systematic errors. 
We summarize our results in section IVll The appendix 
is devoted to a review of dispersive bound. 
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II. BASIC IDEA 

The matrix element of the heavy-to-light semileptonic 
decay B — > nlv is parameterized as 



(ir + (k)\V»(0)\B ( P )) = (p + k 



m r - to; 



/ + (? 2 ) 



(1) 



where — qj^b and q 2 ranges from q^ in — mf(^ 0) 
GeV 2 to q 2 nax = (m^ — mi) 2 . The differential decay rate 
is written as 



this case, the upper and the lower bounds are the solu- 
tions of the quadratic equations whose coefficients are 
determined by the lattice inputs as well as other in- 
puts such as OPE results. For a set of lattice data 
/ = {/+,/°} where /+ = {.f+(q 2 ), ■ ■ • /+(<&)} and 
f° = ' ' ' /°(9jv)}: the upper and lower bounds 

are uniquely determined as a function of q 2 as long as the 
solution of the quadratic equation exists. Let us denote 
the bounds as 

K P /lo(l 2 ->f)> F u P /lo(<l 2 ->f)- ( 3 ) 

Fig. n shows a typical shape of the dispersive bound 



dT(B -> irlv) 
dq~ 2 



24tt 3 



\V uh \ 2 \{vk) 2 ~mlfl 2 \t{q 2 t 



(2) 

Since the discretization error becomes uncontrollable 
for the momentum much larger than Aqcd, the calcu- 
lation of the form factors is feasible only in the large 
q 2 (small recoil) region (q 2 > 16 GeV 2 ), where the spa- 
tial momentum of pion is lower than roughly 1 GeV. The 
CKM matrix clement | V u b | can be obtained by combining 
the experimental data integrated above some q 2 value, 
and the lattice results for the form factor |/ + (g 2 )| 2 in- 
tegrated in the same region with an appropriate kine- 
matical factor. However, due to the limited statistics 
the experimental measurement for this energy range still 
has a large error. Although in principle there is nothing 
wrong with this method, one disadvantage is that the ex- 
perimental data for the rest of the kinematic range with 
much better statistics cannot be used. 

An alternative way is to extrapolate the form factor to 
all kinematic range and combine with the full experimen- 
tal data. For this purpose, we need to have the following 
three ingredients 

1. the experimental data of the partial decay rate 

2 

J q 2 +1 dq 2 dT/dq 2 for the kinematic ranges q 2 < q 2 < 
a 2 

2. the lattice results on the form factor f + (q 2 ) for 
large q 2 region, 



3. a reliable method to extrapolate the form factor to 
region. 



lower q 2 



In order to avoid model dependences, we exploit the dis- 
persive bound to extrapolate the form factor. The disper- 
sive bound is an exact bound on the form factors f (£~ ) 

and f + {q 2 ) for all kinematic range of q 2 II H Ea. 

This is derived from the dispersion relation and the op- 
erator product expansion (OPE) of the two-point corre- 
lation function of the heavy-light (bu) current in deep 
Euclidean region. If we have an additional informa- 
tion on the form factors /°,/ + at some physical kine- 
matic points q\ , • • ■ , q 2 N (from the lattice calculation, 
for instance), the bound can be improved further. In 











- f bound 

- f + bound 

• f + (q 2 ) inputs (no errors) 

• f°(q 2 ) inputs (no errors) 





-25 -15 -5 




15 25 



f° bound 
f + bound 

f (q ) inputs (no errors) 
f° (-q 2 ) inputs (no errors) 




-25 -15-5 5 

q 2 (GeV 2 ) 



15 25 



FIG. 1: bound of f + (t), f°(—t) for one example set of inputs 
/ 

The top figure shows the independent bounds for F^ p n a and 



-F ,, and the bottom figure shows the bounds with the 

up/to ° 

kinematical constraint / + (0) 



/°(0) 
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using a mock data and it is suggested that the inputs at 
three points could give rather stringent constraints. 

One problem is that in practice the additional inputs 
from the lattice calculation always has some theoretical 
errors so that this bound has uncertainties. Lellouch pro- 
posed a statistical treatment and derived a probability 
distribution the dispersive bound 0- He made a ran- 
dom sample set of form factors which obey the Gaussian 
distributions where the central value and the deviations 
are taken from the central values and the error of the lat- 
tice data. For each sample set, the quadratic equations 
which determine the upper and lower bounds are solved. 
He imposed the following condition, 

Condition A (consistency condition): 

• The quadratic equations to determine the up- 
per/lower bounds should have real solutions. 

• The solutions of the upper/lower bounds 
should allow the kinematical condition 
/+(0) = /». 

The sample is accepted when condition A is satisfied and 
discarded otherwise. Then he makes a distributions of 
the upper and lower bounds from only the accepted so- 
lution, which are conditional distributions. 
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FIG. 2: The confidence level bound for /+(*), /°(-t). The 
input data is lattice data(Lellouch) . 

As shown in Fig. the confidence level bound at 
q 2 ~ and q 2 = q^ax nas a huge spread. Although the 
dispersive bound method gives bounds at certain confi- 
dence levels, they are not strong enough to constrain the 
CKM matrix element \V u b\ with high accuracy. 

In the following we would like to consider how to make 
the best use of the experimental data and the theoretical 
results. The large width of the dispersive bound arises 
because the statistical bound was obtained at each q 2 
without considering the correlation in q 2 . However, as 
shown in Fig. 1, each sample of the dispersive bound 
from the Monte Carlo sample based on the lattice data 



gives a very stringent bound and the q 2 dependence can 
be obtained with, say, 10% accuracy. On the other hand, 
from CLEO experiment we also know the q 2 dependence 
up to an overall normalization from unknown factor \V u b\. 
The knowledge of the q 2 dependence of the form factors 
from both the theory and experiment can be an addi- 
tional input for the determination of \V u b\- To extract 
\V u b\, instead of combining the experimental data and 
the lattice results for a single q 2 bin, we make a kind of 
simultaneous fit for all the q 2 bins. Precisely speaking, 
since the lattice results are transformed into the distribu- 
tions of the upper and lower bounds, a little modifications 
are necessary, namely rather than fitting \V u b\ we obtain 
a conditional distributions of | V^e, | using all information 
from the experimental and the lattice results. In other 
words, we can reweight the statistical distribution of the 
form factor bounds using the q 2 distribution from the 
experiment as an additional input, which we will explain 
more in detail in the next section. 

We also note that the soft pion theorem can also be 
used for additional inputs. The soft pion predicts that 
the form factor near the zero pion recoil limit {q 2 — » 
Imax)-! the form factor behaves as 

- t- ( 4 ) 
5 [q,nax) U l-?/m%.' 1 ' 

where g is the B*Bir coupling and f B ~ f B . are the decay 
constants of B and B* mesons [l2|| . Therefore provided 
that g, fs, Jb> are known, the soft pion relation can be 
used to gives additional inputs of the form factor at q 2 nax 
for the dispersive bounds. 



III. METHOD 

In the previous section, we explained the basic idea for 
extracting the CKM element \V u b\ by combing the lattice 
results, experimental data, and the dispersive bound. In 
this section we explain our method more in detail. 

Due to limited statistics and limited energy resolution 
in q 2 etc, usually the whole kinematic range is divided 
into N bin bins (q 2 < q 2 < q 2 +1 ,i = 1, • • • , N bin ) and we 
only know the partial decay rates for those q 2 bins 

r q i+i fir a 2 \v J 2 r q i 

Ff P = / dq 2 ^ = ±J^. / d q 2 \ I+{q 2 )\ 2 X{q 2 f 
J q 2 dq 1 \%2ir A m B J q 2 

(6) 

This partial decay rates are indeed measured by CLEO 
collaboration 

By integrating the dispersive bounds of the form factor 
F^pn (<l 2 ] f) > we can predict the upper and lower bounds 
of the partial decay rate over certain q 2 bins up to the 
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overall factor iVu&l 2 - 



up/lo 

Yi 
G 2 



(/) 
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dt 



where the theoretical upper/lower bounds of the partial 

decay rates are T^ lo (f) = \V ub \ 2 x ^"(f) ■ 

Let us now explain our methods, which is composed of 
five steps. 

1. We generate samples of / from the Gaussian dis- 
tribution with central values and errors from the 
lattice data. The samples of B*Bn coupling 
g and fg ~ /s* are also generated from the 
some distributions based on the best knowledge 
of g and fg, from which the samples of f max = 
{/ + (9Lx)./°(?L)} are produced. They can be 
used as additional data of the form factor for the 
dispersive bounds. 

2. From each sample of {fif ma x\i the dis- 
persive bounds F+p/lo(l 2 ''f>fmax) ( and 

F u P /io(<l 2 'f>fmax)) are derived and each sample 
is accepted or rejected under the Condition A, so 
that a conditional distribution is produced. 

3. For each accepted sample, we compute the up- 
per and lower bounds of the partial decay rate 
l7 ,l °ifi fmaxYs for a given set of q 2 bins. 

4. We create samples of \V u b\ and r^ xp 's. We assume 
\V u b\ distributes uniformly within a conservatively 
wide range, i.e. \V u b\ = [1 x 1CT 3 ,6 x 1CT 3 ]. and 
T^ xp 's are distributed by a Gaussian around their 
central values with variances given by the CLEO 
data. 

5. We further impose the following physical condition 



for the set of samples {/, g, f B , \V ub \, T^} 



CondA- 



Condition B (physical condition): 

• The experimental data r^ p 's should lie 
within the upper and lower bounds from 
the theory simultaneously for all i , i.e. 



< 



\Vub\H° 

(i=l, ••• ,N Un ). 



r 



exp 



< 



\Vub\H 



2 up 



Try this test for all samples, select only 
those combinations which satisfies the condi- 
tion B and reject others. This makes another 
conditional distribution of the set of samples 
{f,9,fB, \Vub\,Tl xp } C ondAB which satisfies both 
Condition A and B. 

Mathematically the original probability for the set of 
samples {/, g, fg, \V ub \, T exp } is given by the following 
product 

P(f,gJ B ,\Vu b \,T exp ) 
= NV f (f)V g (g)V fB (fB)Ve Xp (r exp )V lVubl (\V u b\),(8) 



where N is the normalization constant. The probability 
V j!,Vf B ,P exp are Gaussian distributions based on theo- 
retical or experimental data, where as Pg,P\v ub \ are uni- 
form distributions. 

After imposing the Condition A, the probability be- 
comes 

P A (f,9,f B AV u blT e * p ) 

_ J N'P(f,g,fg,\V u b\,T exp ) (if CondA is satisfied), 
(otherwise), 

where N' is another normalization constant to make the 
total probability unity. By further imposing the Condi- 
tion B, the probability then becomes 

P A B(f,9jB,\Vub\,T eXp ) 

N"P(f,g,f B> \V u b\,T exp ) (if Cond A+B is satisfied), 
(otherwise), 

where N" is another normalization constant. 

Using this last conditional distribution Pab , the distri- 
butions of any quantities defined from the set of samples 
{/; <?> /si \ V u b\, F exp }AB, can be derived, which of course 
includes that of \V ub \- 



IV. RESULTS 

A. The setup 

In this section we explain our results based on the lat- 
tice results and experimental data. The experimental 
data are taken from the results from CLEO collabora- 
tion Q • Table |U shows the partial branching fraction 
Fj/r totQ ;. In order to convert Ti we take r tota / = [4.29 ± 



g 2 bin 


[(/ dr/dq 2 )*dq 2 ]/T total 


< q 2 < 8GeV 2 
8GeV 2 < q 2 < l&GeV 2 
WGeV < q 2 < q 2 max 


(0.43 ±0.11) x 10" 4 
(0.65 ±0.11) x 10" 4 
(0.25 ± 0.09) x 10" 4 



TABLE I: CLEO data of partial branching ratio for B — > ttIu 



0.04] x 10~ 13 GeU 0. There are lattice calculations from 
four different lattice groups in the quenched approxima- 
tion 0, 0, |23, 0| • In this section, we take the lattice 
result from JLQCD collaboration [T3I - We picked the 
form factor data at three points of q 2, s for our analysis, 
which are listed in Table ITT1 The recent B*Bir coupling 
g in quenched lattice QCD gives g = 0.48 ± 0.03 ± 0.11 
for the static limit and g = 0.58 ± 0.06 ± 0.10 for the 
b quark from the quenched lattice calculations 13] . The 
QCD sum rule gives g = 0.36±0.10 Q. The predictions 
from quark models are in the range 0.3 < g < 0.8 0]- in 
our analysis we took uniform distribution g = [0.3,0.9]. 
The decay constant fg is known with 30% error. Since 
this error is much larger than the splitting of fg* and fg, 
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q 2 (GeV 2 ) 


f + (l 2 ) 


/V) 


17.79 
19.30 
20.82 


1.03 ± 0.22 
1.24 ±0.21 
1.54 ±0.27 


0.407 ± 0.092 
0.45 ±0.11 
0.51 ±0.14 



TABLE II: lattice results of f + {q 2 ) and f°{q 2 ) by JLQCD 



we took /s» = Jb = 190 ± 30 MeV for simplicity, taking 
the best estimate by the lattice calculation UjJ . 

In step 1 and 4, based on the above information we 
generated O(10 7 ) set of samples for {/, g, fg} , 2,000 set 
of samples for {T exp }, and 2,000 samples for \V u b\. 

Imposing the condition A in step 2, 2000 set of sam- 
ples of {/, g, fs} survives. In Fig. El the histogram for 
the samples of / is shown. We find that the dispersive 
bound actually reduces the variance of the lattice data 
significantly. 



noted the open window for each sample from F^ p , lo (q 2 ) is 
very small. Integrating the squares of these upper/lower 
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FIG. 3: distribution of form factor f+(q 2 = 17.79GeF 2 ). The 
Top figure is the original samples from Gaussian distribution 
and the bottom is the samples which survived condition A. 

In step 3, we compute the upper/lower bounds for the 
partial decay rate divided by \V u b\ 2 for each q 2 bins. To 
illustrate how it works, we pick up 3 generic set of samples 
( samples 1, 2 and 3 ). In Fig.EJwe show the upper/lower 
bounds of the form factors for each sample. It should be 



FIG. 4: Upper/lower bounds of the form factors / + for sample 
1,2,3. 
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FIG. 5: Upper/lower bounds of the partial branching frac- 
tion over the three energy bins for sample 1, 2, 3, where 
gup/la _ j^ p ^ lo ^v u b\ 2 /Ftotai- While the overall normaliza- 
tion is unknown up to |Vui>|, we have set \V u b\ = 3.71 x 10~ 3 
for illustration. 

bounds for the form factors for each sample, we obtain 
the upper/lower bounds for ^ p ^ l ° and the partial branch- 
ing fraction B™ p/l ° = ^ p/l °\V ub \ 2 /T tota i for each q 2 bin 
as shown in Fig. [SJ For comparison we also show the 
the partial branching fraction by CLEO collaboration in 
Fig. El 

In step 5, we impose condition B. As can be seen from 
Fig. Eland Fig. El it is expected that the conditional prob- 
ability with condition A±B would be highly suppressed 
for samples 1 and 3 but not so suppressed for sample 2. 
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0.8 - 




FIG. 6: The partial branching fraction Bi's over the three 
energy bins from CLEO, where Bi = Yi/Y to tai- 

The histogram for the samples of / after imposing con- 
dition B is shown in Fig. We find that the condition B 
reduces the variance of the form factor distribution even 
further. In order to understand the mechanism of the re- 



FIG. 8: scatter plot of 73/72 vs 72 and 71/72 vs 72 
under condition A. As a reference we also show the favored 
region of Yf p /T^ v and r^ p /r^ p 
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FIG. 7: distribution of /+ (q 2 = 17.79GeV 2 ) from the samples 
which survived condition A and B . 



duction of the variance, we show show the scatter plot of 
73/72 vs 72 or 71/72 vs 72 in Fig. El - We find that there 
are clear correlations. From the CLEO data in Fig. El 
it can be shown that the experimental data favors the 
range T e 3 xp /T e 2 xp ~ 0.2 - 0.5 and T exp /T e 2 x p ~ 0.4 - 0.9. 
Since the open window of upper/lower bounds for each 
sample is quite small, imposing condition B is similar to 
restricting 73/72 = 0.2 — 0.5 and 71/72 = 0.4 — 0.9 in 
Fig. |H1 Since there is a strong correlation this restricts 
72. In fact after imposing condition B, the scatter plot 
corresponding to Fig. El looks as Fig. El 



We give our final results for various physical quanti- 
ties by looking their distributions which can be extracted 
from the conditional distribution Pab- 

First, the dispersive bound is shown in Fig.^3 Fig. 1111 
shows the histogram of \V u b\. From this distribution we 
obtain 

\V ub \ = [3.73 ±0.53] x 10~ 3 . (9) 

This can be compared to the \V u b\ determination using 
only the raw lattice data and the CLEO result at highest 
q 2 bin. For example, if we use the JLQCD results of 
the differential decay rate for q 2 > 16 GeV 2 , we obtain 
\V ub \ = [3.1 ±0.9] x 10~ 3 . This means that the \V ub \ error 
of 30% in the conventional method is reduced to 14 % in 
our new method. 

This error reduction is rather remarkable and requires 
some explanation on which part of the analysis con- 
tributed most. One important ingredient is the condition 
A and the other important ingredient is the condition B. 
Also, whether or not we have a soft pion input is another 
important point. In order to see the effect of each ingre- 
dients, in addition to our full analysis we also carried out 
analyses with the following three cases. 

1. Condition A only without soft pion input, 

2. Condition A only and with soft pion input, 

3. Condition A + B but without soft pion input. 
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FIG. 9: scatter plot of 7 f /j 2 h vs and it lit vs 7 | h 
after imposing condition A+B. As a reference we also show 
the favored region of Y^ v and V^ v jY^ v . 
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2 

FIG. 10: CLB for f + (q 2 ). Here we used JLQCD's lattice 
input, and CLEO's experimental data. 



Fig. El Fig. El and Fig. El show the corresponding 
results. 

The figures suggests that both condition A and condi- 
tion B reduces the error. In fact, the corresponding \V u b\ 
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FIG. 11: Histograms of \V u b\- The original distribution was 
assumed to be flat in the range |Ki>| = [1 5 6] -3 (left). After 
imposing condition B, the resulting distribution distributes 
with small variance. 

values are 

\V ub \ = [2.92±°;<*] x 1(T 3 for case 1, (10) 

1Kb | = [2.99±°-j|] x 1(T 3 for case 2, (11) 

1Kb | = [3.61±°;g] x 1(T 3 for case 3, (12) 

so that the error for \V ub \ is has 20%, 20% and 14% fra- 
cases 1, 2 and 3 respectively. Thus we find both condition 
A and condition B reduces the error and having both of 
them together gives the significant reduction. It is also 
found that the input from the soft pion theorem does not 
change the result of |K&| so much. 

We also obtain product we also obtain bounds at 66% 
confidence level (66% CL) for / + (0) as a byproduct 

0.126 < /(0) < 0.293. (13) 

which will be very useful for predicting the two body 
decay rate in QCD factorization . 

V. SYSTEMATIC ERRORS 

In this section we discuss possible systematic errors in 
our analysis. 
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FIG. 12: the case of Condition A only, without soft pion 
input. 
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FIG. 13: the case of Condition A only, with soft pion input. 



The main systematic error arise from the lattice data 
and from the experimental data. These errors are of 
course already included in the errors of our input data. 
However, it would be nice if we could estimate these er- 
rors in a different fashion. One way to do so is to use the 
input data from other groups which has different sys- 
tematic errors. In the following we present our result 
using different lattice inputs. Although we should also 
perform analyses using other experimental data such as 
those from Belle or BaBar, this will be left as a future 
problem since there is no publicly available data. There- 
fore we will not discuss the systematic errors from exper- 
imental data in this paper. 

There are four quenched lattice calculations with dif- 
ferent systematic errors. JLQCD collaboration used non- 
relativistic effective theory for the heavy quark on a 
coarse lattice, in which the dominant systematic errors 
are the momentum dependent discretization error and 




/ ! 
; i 

i 



-25.0 -15.0 -5.0 5.0 15.0 25.0 

FIG. 14: the case of Condition A+B, without soft pion input. 

chiral extrapolation error. Fermilab collaboration also 
used a different type of nonrelativistic effective theory 
on relatively coarse lattices, using different chiral extrap- 
olation method. APE collaboration and UKQCD collab- 
oration used relativistic action for the heavy quark on 
finer lattices, in which the dominant error is the presum- 
ably momentum independent discretization error from 
the large heavy quark mass in lattice units. In Table ITTT1 
we give the input parameters from APE, UKQCD and 
Fermilab collaborations. The setups for the numerical 



collaboration 


q 2 (GeV 2 ) 


f + (l 2 ) 


/V) 


Fermilab 


17.23 
20.35 
23.41 


1 -in + 0.29 

r.io_ 19 
-, 79 +0.31 
L - ' z -0.27 

2 10+ ' 43 


0.83i° ;£ 
1.001°;?° 


APE 


13.6 
17.9 
22.1 


70+° 13 

u - ' u -0.09 

1 n5+ 015 

1 qfi+°- 24 
r.yo_ 30 


0.46l°o°i 

o.59tS:SI 

0.80l°; 06 3 


UKQCD 


16.7 
18.1 
22.3 


0.90 ±0.22 

1 io+°- 28 

1 - iu -0.22 


0.57l°;°« 
0.611°-?? 
0.79+°-?? 



TABLE III: lattice results of f + (q 2 ) and f°(q 2 ) by Fermilab, 
APE and UKQCD. 

analysis are almost the same except for the number of 
generated samples for {f,g, Jb} in step 1. In the anal- 
yses using Fermilab, APE and UKQCD data, we have 
generated 0(1O 5 ) samples. The results using Fermilab 
data are 

1Kb | = [3.59±g;g] x 1(T 3 , (14) 
0.128 < /(0) < 0.339 (at 66% CL), (15) 

those using APE data are 

I Kb | = [3.35t{j:2jj] x 10- 3 , (16) 
0.132 < /(0) < 0.302 (at 66% CL), (17) 



9 



those using UKQCD data are 



1Kb = [3.37+^] x 10~ 6 , 
0.154 < /(0) < 0.353 (at 66% CL). 



(18) 
(19) 



Another possible source of systematic errors are the 
theoretical input values of Xz-(O), Xt(0) for the disper- 
sive bound which is obtained by the operator product 
expansion (OPE) using perturbative QCD and some es- 
timate of vacuum condensation values. (See appendix 
for more details.) As is explained in the appendix these 
input values are 



XT = [5.60 ±0.17] x 10~ 4 , 
XL = [1.50 ±0.03] x 10~ 2 , 



(20) 
(21) 



where the errors are from unknown 2-loop perturbative 
corrections and uncertainties in the vacuum condensa- 
tion values. By changing the input values of xt and xl 
by 1 er, the central values of |Kb| changes by 2% when 
JLQCD data are used, which is negligible compared to 
other errors. We also found that the same is true when 
Fermilab, APE, or UKQCD data are used. 

Yet another systematic error is the quenching error in 
the lattice results. Since there is no way to control this 
error other than performing unquenched lattice calcula- 
tions, which is beyond the scope of this paper. We leave 
it an open question and wait for the unquenched lattice 
results, which is expected to appear near future. 



VI. CONCLUSION 

In this paper, we have proposed a method to determine 
\V u b\ by combining lattice results, dispersive bounds, and 
experimental data. Based on Lellouch's idea we consid- 
ered the statistical distributions of the form factors and 
their bounds. Our proposal is to restrict the distribution 
by imposing a physical condition using the experimental 
data. This gives a significant reduction in the errors of 
\V u b\ ■ As as result we obtained 



IK, 



[3.73 ±0.53] x 10" 



(22) 



As a by product we also obtained a 66% confidence level 
bound for the form factor at small q 2 , 



0.293 < /(0) < 0.126 



(23) 



which will be very useful for predicting the two body 
decay rate in QCD factorization. 

The recent Belle results of |Kfe| from B — > X u lv [TT|is 



IK 



[4.66±0.28±0.35±0.17±0.08±0.58] x 10~ 3 , (24) 



where the errors is the statistical, systematic, b — > c, 
b — > u, and theoretical. The lattice QCD based |K&| 
values by CLEO collaboration from B — > irlv is 



IK 



[2.88 ± 0.55 ± 0.30tg;fJ ± °- 18 ] x 10" 3 ,(25) 



where the errors are statistical, experimental systematic, 
theoretical and plv form factor shape. It should be noted 
that this discrepancies are reduced by our method al- 
though the input lattice data are the same. 

For future studies, we should elaborate our method in 
such a way that we treat the systematic and statistical 
errors separately which would give a more careful anal- 
ysis of errors. This is in principle straightforward but 
requires more numerical calculation since the generation 
of random samples should be much increased. 

Another improvement is to update the input data of 
form factors from lattice, especially the unquenched one, 
and the input data of partial decay rates from B factories 
such as Belle and BaBar. Both of these will be available 
near future. 



Appendix: Review of the Dispersive Bound 

In this appendix, we briefly summarize the disper- 
sive bounds on the form factors, which is fully discussed 
in Ref. [2|. The matrix element of the heavy-to-light 
semileptonic decay B — > irlv is parameterized as Eq. . 
The form factors / + and /° satisfies the following kine- 
matical constraint 



/+(0) = /». 



(26) 



In order to derive bound on f + {q 2 ) and f°{q 2 ) we con- 
sider the following two-point function 

n^(q) = i J d 4 xe iqx (0|TV' i (x)V l/ (0)|0) 

= -l^V - q^l^Tiq 2 ) + -f -f n L (g 2 )(27) 

where Ht(l) corresponds to the contribution from the 
intermediate states with J p = 1~(0 + ). In the deep Eu- 
clidean region which is far from the physical cut, this two 
point function can be evaluated by perturbative QCD re- 
liably. 

Inserting all possible intermediate states \T) we obtain 
the following result for the imaginary part of : 

-(g^q 2 - g"g")Imn T + <fq v Imn. L 

= i^(2 7 r)^ 4 ( (Z -p r )(0|^|r)(r|K t |0), (28) 
r 

We now restrict ourselves to include only contributions 
of the Bit and the B* states. Since each hadron state 
contributes positively to the spectral function, we derive 
the following two bounds for the form factors 



Of f 

Imn L (u > ^L^y/(t-t+)(t-t-) 



32tt 



t 3 



9(t-t+), 



ImII T (t) > 7T 



5{t-m B , 2 ) 



1 [(t-t + ){t-t_)f/ 2 



32tt 



f + (t)\ 2 e(t~t + ), 
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where t = q 2 and t± = (ms± m^) 2 ■ In the above 
equations we have limited ourselves to contribution of 
the B* state with 



(0\V\B*(r,p)) 



h 



(29) 



In QCD, the polarization functions II ^ and Tlx obey once 
and twice subtracted dispersion relation(q 2 = — Q 2 ), i.e.: 



Xl(Q 2 



xt(Q 2 



d 



(-Q 2 u L (Q 2 )) 



dQ 2 ' 

dttlmn L (t) 



o *{t + Q 2 ) 2 



1 / d 



2 V dQ 2 



(-Q 2 U T {Q 2 )) 



dt tImYl T (t) 

V(t+Q 2 r 



(30) 



(31) 



Thus it is now possible to get bounds on the form factors 
to inserts Eas. l|30l31jl to Eas. l|29l29|) 



Jl(Q 2 ) > - / dtk L (t,Q 2 )\f (t)\ 2 (32) 
Jt(Q 2 ) > - / dtk T (t,Q 2 )\f + (t)\ 2 (33) 



where J£(t))&l(t) are 

Jl = xl(Q 2 ), 
Jt = Xt(Q 2 ) 



k L = 



(m B . + Q 2 ) 3 V fa 

3 [(t-t + )(t-t_)] 1/2 M- 
32^ (* + Q 2 ) 2 t 2 

3 [(t-t + )(t-t_)] 3/2 
96tt (t + Q 2 )t 2 



(34) 

' * , (35) 
(36) 
(37) 



To obtain the bounds of the form factor for values of t 
in the range [0, t_], we map the complex t-plane into the 
unit disc in complex z-plane with the conformal transfor- 
mation 



l + z _ (t+-t)(t--t) 



1 - z 



U - t- 



This maps the contour along the physical cut in 
Eas. 1)32133(1 onto a unit circle in z-plane so that 



Jl(Q 2 ) > 
MQ 2 ) > 



dz 
2niz 

dz 
2niz 



\Mz,Q 2 )fo(z)\ 2 , (38) 



b T {z,Q 2 )f+{z)\ 2 , (39) 



where we have used the fact that k(t, Q 2 ) is positive def- 



inite quantity. The functions <pL{T) are defined as 



t+ l + z 



t+ l + z 



3t + t_ 1 l + z 
4tt t + - tl (1 - zfl 2 \yt + -t- ' 1-: 

, = / 1 (1+^) 2 / 

9T \ TT{t+-t_){l- Zf/ 2 \y t+-t_ ' 1-Z 

( j t + -Q 2 i+ z y 3 
{it^r + T—z) < 

and satisfy cf>(z, Q 2 ) > for t 6 [0, t~]. For simplicity, we 
define an inner product on the unit circle 



(g\h) = 



dz 



z\ = l 2tTIZ 



g*(z)h(z), 



(40) 



so that the inequality Eq. I|38|l . Eq. (|39|l can be written 

J > (0/10/)- (41) 
And define the function gt(z) and the matrix M(f(i)) 



as 



9t{z) = 



1 - z*(t)z' 



MWt)) - »» fM) 
\(9tW) (9t\9t),J 

then if f(z(t)) has no pole in the range [0, t~], f{z(t)) 
can be related the inner-product of gt and <j)f as 



iatW) =<p(z(t),Q 2 )f(z(t)) 

Now, because positivity of inner products 
detM(/(f),/(ti)) >0 



(42) 



(43) 



is satisfied, by eliminating (<fcf\4>f) with I|38I39|I we get 
the form factors bound as 



l/WI 2 < J(Q 2 ) 



l 



i 



l-\z(t)\ 2 cf>(z(t),Q 2 ) 



(44) 



If we want to get more conditionally bound, using some 
values of the form factors at t\,t2, • • • , th and defining 
the matrix M(f(t), f) as 



M(f(t),f) 



(WW) W\gt) (Mint*) ■■■ W\9t N )\ 
(9t\4>f) (gt\gt) (gt\gt 3 ) ■■■ (gt\g tN ) 
(9t 3 \(t>f) (gt 3 \gt) (gt 3 \gt 3 ) ■■■ (gt 3 \g tN 



\(9t N \^f) (gt N \gt) (gt N \g t3 



(9t N \gt N )J 
(45) 



11 



where / stand for {f(z(t 1 ))J(z(t 2 )),---J(z(t L ))}. 
This matrix M(/(t), /) satisfies 

detM(f(t),f)>0, 



By eliminating (4>f\4>f) using (|38I39[1 . the inequality 146|) 
leads to bounds on the form factors f(t), 



FUt;f)<f(t)<F up (t;f). 



This upper and lower functions F, 



up{lo) 



stand for 



F, 



up/lo 



-M + / - vW) ■ A(t) 

a ■ <j){t, Q 2 ) 



(46) 



where a, f3(t), and ^c(Q 2 ) ■ A(t) are known functions of 
Q 2 , ti, and f(z(ti)) (i = 1, • • • ,L). The particular form 
of the Eqs. P^l can be found in Ref. 0- 

Now, if f+(t) has a single pole at t — t p away from the 
cut (in fact t p = ms* 6 [£_,£+]), Eqs. (|42|l becomes 



So we can make kinematical constraint dispersive bound 
like bottoms of There may exist 4 pattern no kine- 
matical constraint bound, from there kinematical con- 
straint bound can be made respectively same method. 

In this calculation, they use OPE and calculate the 
coefficient by perturbative QCD. In this paper, we take 
Xt(0) = (5.60±0.17)x1CT 4 ,xl(0) = (1.50±0.03) x 1CT 2 
(See Ref. |2j). The values of \ in MS scheme are given in 
Table. IIVI These coefficients are obtained from the op- 
erator product expansion of hadron current correlators 
in perturbative approximation. In order to control the 
higher order contributions which grow for larger — Q 2 , 
we take Q 2 = 0. The table shows that the 0(a s )/lLoop 
are 17%, 7% for xt,Xl so that the perturbative expan- 
sion is under control. The unknown 2-loop errors can be 
estimated by either naive order counting or by compar- 
ing results with different renormalization schemes, which 
amount to less than or equal to 1%. We therefore esti- 
mate the uncertainties in \ from the condensates, which 
are of order 2% and 3% for xl and xt- 



Res((j)f;z(tp)) 
z(t p ) - z(t) 



so {gt\<t>f) have the single pole at z = z(t p ). In this case 
we define the function 4> p (z), which satisfies \<p p \ = at 
\z\ = 1, as, 



z{t p ) 



l-z*(t p )z 



This function <j) p f p dose not have a pole at z — z(t p ), so 



((f>pfp\9t) 
(<l>pfp\<l>pfp) 



(4>f\4>!) < J(Q 2 ) 



are satisfied. Thus, the bound of Eqs. 146|) holds 

*?„(*;/)</(*) <*£,(*;/). 

where Ff up are the functions Fi Q>up obtained by replac- 



ing (j> with 4> p . In this paper we denote 



4>i Fi 0iUp for /+. In this way, we get the the independent 
bounds on the form factors using the given values of the 
discrete set of points ,for example top figure of Fig. 
with no kinematical constraint. But form factors not in- 
dependent but satisfy the kinematical constraint Eas l26l 
So we make a new dispersive bound, use 



f + (0)u P = F° (0,/) 



(47) 



as the additional input of the new bound for f + (t), 
and consider this upper bound as new appear bound for 
f + (t). And equally make a new lower bound for f + (t), 
using the additional input 



f + (0) lo = F+(0,f). 



(48) 





Xl(Q 2 = 0) 


Xt(Q 2 = 0) 


lLoop 


1.3 x 10" 2 


5.1 x 10 -4 


0(a s ) /lLoop 


17% 


7% 


total pert 


1.5 x 10" 2 


5.5 x 10 -4 


m u < uu > /lLoop 


-2% 


3% 


< a s G 2 > /lLoop 


0.04% 


-0.05% 


total 


1.5 x 10" 2 


5.6 x 10~ 4 



TABLE IV: perturbative and non-perturbative contributions 
to subtracted polarization function xl and xt by MS- 
scheme. "1 Loop" stands for the 1-loop result and "totalpart" 
stands for the results through 0(a s ). m u < uu > /lloop,< 
G 2 > are the higher twist contributions from the conden- 
sates. 
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